egrandis = read.table("egrandis.csv", head=T, sep=";") ##Lendo a tabela egrandis.csv

par(mfrow=c(2,1)) #Ajuste para que os prximos grficos sejam mostrados um por pgina

boxplot(dap~regiao, data=egrandis) 

boxplot(dap~rotacao, data=egrandis) 

par(mfrow=c(1,1)) #Ajuste para que os prximos grficos sejam mostrados um por pgina 


#Normalidade da altura do conjunto total de rvores com um grfico quantil-quantil contra a distribuio normal.
qqnorm(egrandis$ht) 
qqline(egrandis$ht) 
